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This  paper  develops  a  three-dimensional  analytical  model  of  a  fluid-loaded  acoustic 
coating  affixed  to  a  rib-stiffened  plate.  The  system  is  loaded  by  a  plane  wave  that  is 
harmonic  both  spatially  and  temporally.  The  model  begins  with  Navier-Cauchy 
equations  of  motion  for  an  elastic  solid,  which  produces  displacement  fields  that  have 
unknown  wave  propagation  coefficients.  These  are  inserted  into  stress  equations  at  the 
boundaries  of  the  plate  and  the  acoustic  coating.  These  stress  fields  are  coupled  to  the 
fluid  field  and  the  rib  stiffeners  with  force  balances.  Manipulation  of  these  equations 
develops  an  infinite  number  of  indexed  equations  that  are  truncated  and  incorporated 
into  a  global  matrix  equation.  This  global  matrix  equation  can  be  solved  to  determine 
the  wave  propagation  coefficients.  This  produces  analytical  solutions  to  the  systems’ 
displacements,  stresses,  and  scattered  pressure  field.  This  model,  unlike  previously 
developed  analytical  models,  has  elastic  behavior  and  thus  incorporates  higher  order 
wave  motion  that  makes  it  accurate  at  higher  wavenumbers  and  frequencies.  An 
example  problem  is  investigated  for  three  specific  model  results:  (1)  the  dynamic 
response,  (2)  a  sonar  array  embedded  in  the  acoustic  coating,  and  (3)  the  scattered 
pressure  field.  An  expression  for  the  high  frequency  limitation  of  the  model  is  derived. 
It  is  shown  that  the  ribs  can  have  a  significant  impact  on  the  structural  acoustic 
response  of  the  system. 

Published  by  Elsevier  Ltd. 


1.  Introduction 

Underwater  vehicle  hulls  are  typically  metallic  with  acoustic  coatings  that  are  made  of  a  soft  polymer  material.  These 
hulls  are  usually  internally  reinforced  in  one  direction  to  provide  increased  stiffness  against  the  hydrostatic  forces  that  act 
on  the  hull  when  the  vessel  is  submerged.  Acoustic  coatings  are  used  to  produce  quieter  vehicles,  contain  sonar  sensors, 
inhibit  drag,  and  prevent  biological  fouling  on  the  exterior  of  the  vehicle.  Understanding  the  structural  response  of  such  a 
system  is  important  for  the  design  of  new  underwater  vehicles  and  analysis  of  existing  underwater  vehicles.  An  analytical 
model  with  this  type  of  configuration  will  allow  predictions  of  the  dynamic  response  of  an  acoustic  coating  on  a  marine 
structure  and  a  sonar  system  embedded  in  the  coating.  Over  the  years  there  has  been  a  progression  of  modeling  these 
systems,  from  thin  unreinforced  structures  to  thick  unreinforced  structures  and  thin  reinforced  structures.  The  specific 
problem  of  an  analytical  model  of  a  thick  reinforced  structure  has  not  been  previously  addressed. 

The  research  of  plate  theory  has  been  an  ongoing  field  of  study  for  many  years.  Thin  beam  (Euler-Bernoulli)  theory 
originates  in  the  eighteenth  century  [1].  This  model  is  sometimes  called  a  flexural  wave  model,  and  an  extremely  similar 
model  exists  in  plate  theory.  Mindlin  [2]  modified  thin  plate  theory  to  include  the  dynamics  of  rotary  inertia  and  shear 
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effects,  and  this  extended  the  low  frequency  and  wavenumber  region  of  the  model.  Analysis  of  wave  propagation  in  fully 
elastic  plates  consisting  of  one  or  multiple  layers  using  the  Navier-Cauchy  equations  of  motion  has  been  studied 
extensively  and  is  well  documented  [3].  Fluid  loading  has  been  added  to  thin  plate  analysis  [4]  and  thick  plate  single  and 
multilayer  analysis  [5].  Using  these  techniques,  the  dynamic  response  of  fluid-loaded  acoustic  coating  on  an  unreinforced 
plate  can  be  calculated  and  analyzed  for  both  thin  and  thick  structures. 

Over  the  years,  thin  plate  theory  has  been  expanded  to  include  stiffening  effects  of  ribs  that  are  typical  of  many  marine 
and  aviation  structures.  These  ribs  are  almost  inclusively  placed  on  one  side  of  the  plate  to  reinforce  the  structure. 
The  response  of  a  periodically  supported  beam  to  a  load  with  fixed  wavenumber  and  frequency  has  been  investigated  [6]. 
The  response  of  periodically  stiffened  fluid-loaded  plates  to  harmonic  loading  [7]  and  line  and  point  forces  [8]  has  been 
established.  The  problem  of  aperiodicity  in  the  stiffeners  has  been  solved  [9].  This  problem  was  also  investigated  for  a  finite 
number  of  equally  spaced  stiffeners  [10]  and  randomly  spaced  stiffeners  [11].  Asymptotic  models  of  plate  radiation  into 
fluid  fields  have  also  been  developed  [12].  The  problem  of  an  arbitrary  number  of  dissimilar  ribs  which  are  randomly 
spaced  has  been  solved  [13].  A  model  of  radiation  of  sound  from  a  reinforced  plate  subjected  to  a  line  load  has  been 
developed  [14].  In  this  paper,  the  author  incorporated  the  size  of  the  ribs  into  the  analysis  of  the  radiated  sound  field. 
Sound  radiation  from  a  plate  reinforced  by  infinite  sets  of  orthogonal  ribs  has  been  studied  [15].  A  similar  problem  with 
doubly  periodic  orthogonal  ribs  has  been  studied  with  respect  to  the  diffraction  of  an  incoming  plane  wave  [16].  It  is  noted 
that  these  papers  with  reinforcing  ribs  [6-16]  use  some  form  of  thin  plate  theory,  and  the  resulting  frequency  limit  where 
the  model  assumptions  are  valid  depends  on  the  thickness  of  the  plate,  but  for  most  plates  it  is  typically  in  the  hundreds  of 
Hertz.  These  thin  plate  theories  are  based  on  a  single  flexural  wave  traveling  in  the  plate,  and  do  not  include  the  dynamics 
of  higher  order  wave  motion. 

The  acoustic  response  of  a  fully  elastic  cylindrical  shell  with  a  complete  acoustic  coating  has  been  researched  [17]  as 
well  as  that  of  a  cylindrical  shell  with  a  partial  acoustic  coating  [18].  Rib  effects  are  not  present  in  either  of  these  papers. 
Higher  frequency  analysis  is  possible  using  numerical  methods  such  as  finite  element  analysis  [19];  however,  these 
computations  can  be  time  consuming  and  frequently  have  stability  problems  when  Poisson’s  ratio  of  the  coating  material 
approaches  one  half.  Recently,  a  fully  elastic  solution  to  a  thick  plate  containing  discrete  masses  was  developed  [20]. 
This  method  can  be  extended  to  model  the  behavior  of  an  acoustic  coating  on  a  ribbed  backing  plate. 

This  paper  derives  an  analytical  model  of  a  reinforced  elastic  plate  system.  This  model  was  specifically  developed  so 
that  this  system  could  undergo  analysis  at  higher  wavenumbers  and  frequencies.  Previously,  only  the  low  wavenumber 
and  low  frequency  region  of  this  system  was  available  for  analytical  modeling.  This  is  of  interest  because  many  sonar 
systems  operate  at  frequencies  well  above  the  limits  of  flexural  wave  dynamics  and  if  the  structure  is  reinforced, 
traditional  thick  plate  theory  does  not  include  the  dynamics  of  the  stiffeners.  If  the  system  has  fairly  rigid  stiffeners,  then 
their  effects  must  be  incorporated  into  the  model.  The  work  here  was  developed  so  that  the  structural  acoustic  analysis  of 
underwater  vehicles  could  be  accomplished  at  elastic  wavenumbers  and  frequencies  using  an  analytical  model. 

The  modeled  system  consists  of  a  plate  in  contact  with  a  fluid-loaded  acoustic  coating  on  one  side  and  a  series  of 
equally  spaced  ribs  on  the  other  side.  The  structure  is  loaded  via  the  fluid  with  an  incoming  acoustic  wave.  The  plate  and 
the  acoustic  coating  are  modeled  as  three-dimensional  fully  elastic  solid  bodies,  the  fluid  is  modeled  as  a  three- 
dimensional  acoustic  field,  and  the  ribs  are  modeled  using  the  Timoshenko  beam  equation,  the  torsional  wave  equation 
and  the  bar  wave  equation.  The  formulation  begins  with  elasticity  theory  where  the  motion  in  the  plate  and  the  acoustic 
coating  are  modeled  as  a  combination  of  dilatational  and  shear  waves.  These  waveforms  can  be  used  to  determine  the 
three-dimensional  displacement  fields  with  unknown  coefficients.  These  displacements  are  inserted  into  stress  and 
continuity  equations  at  the  system  boundaries  and  interfaces  that  contain  the  system  excitation,  the  fluid  loading,  and  the 
force  of  the  ribs  on  the  structure.  Using  an  orthogonalization  procedure  produces  an  m-indexed  mathematical  model  of  the 
system  where  each  m  index  is  a  set  of  equations  12  rows  by  an  infinite  number  of  columns.  All  of  the  m-indexed  equations 
can  be  combined  and  this  yields  a  matrix  system  of  infinite  extent,  which  is  truncated  to  a  finite  number  of  terms.  Inverting 
this  matrix  solves  for  the  unknown  coefficients  and  produces  system  displacements,  stresses,  and  pressure  field  solutions. 
The  elastic  model  is  compared  to  a  previously  developed  thin  plate  model  at  low  frequency  and  wavenumbers  to  insure 
accuracy  and  consistency  to  prior  work.  A  numerical  example  problem  is  then  studied  with  specific  interest  in  dynamic 
response  of  the  structure,  reception  of  acoustic  signals  and  strength  of  the  scattered  field.  The  frequency  limitation  of  the 
new  elastic  model  is  discussed.  It  is  shown  that  the  ribs  can  significantly  affect  the  system  response. 


2.  Governing  equations 

The  system  model  is  that  of  a  fluid-loaded  solid  layer,  called  the  acoustic  coating,  in  contact  with  a  rib-stiffened  solid 
layer,  called  the  backing  plate,  as  shown  in  Fig.  1.  This  problem  is  analytically  modeled  by  assuming  the  pressure  in  the 
fluid  is  governed  by  the  acoustic  wave  equation,  both  solid  layers  are  governed  by  fully  elastic  equations  of  motion,  and  the 
rib  stiffeners  are  modeled  using  the  Timoshenko  beam  equation  for  the  bending  motion,  the  torsional  wave  equation  for 
the  twisting  motion,  and  the  bar  wave  equation  for  the  longitudinal  motion.  The  acoustic  coating  is  excited  on  the  top 
surface  with  an  incoming  acoustic  wave.  The  ribs  on  the  bottom  of  the  backing  plate  are  equally  spaced  at  a  distance  of  L  in 
the  x-direction.  The  model  uses  the  following  assumptions:  (1)  the  excitation  or  forcing  function  acting  on  the  plate  is  a 
plane  wave  at  definite  wavenumbers  in  the  x-  andy-directions  and  frequency  in  time,  (2)  the  displacement  in  both  plates  is 
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three-dimensional,  (3)  both  plates  have  infinite  spatial  extent  in  the  entire  x-  and  y-directions,  (4)  the  fluid  medium  has 
infinite  spatial  extent  in  the  entire  x-  and  y-directions  and  infinite  spatial  extent  in  the  positive  z-direction,  (5)  the  interface 
of  the  acoustic  coating  and  the  backing  plate  have  equivalent  displacements  and  stresses,  (6)  the  motion  of  the  ribs  is 
independent  in  three  directions  and  is  composed  of  twisting  along  the  y-axis,  extension  in  the  y-direction,  and  bending  in 
the  z-direction,  and  (7)  the  particle  motion  is  linear,  and  (8)  the  fluid  medium  is  lossless. 

The  motions  of  both  the  backing  plate  and  the  acoustic  coating  are  governed  the  Navier-Cauchy  equations  of  motion 
written  in  vector  form  as 


n' V2u(x,y,z,  t) + a + jU)VV  •  u  (x,y,z,  t)  =  p  ,  ( 1 ) 

where  p  is  the  density,  X  and  p  are  the  (usually  complex)  Lame  constants,  t  is  time,  •  denotes  a  vector  dot  product, 
and  u(x,y,z,t)  is  the  three-dimensional  Cartesian  coordinate  displacement  vector.  Using  standard  techniques  [21]  to  solve 
Eq.  (1)  and  incorporating  the  periodicity  of  the  system  [6]  into  the  solution  yields  the  displacements  as 

m  =  +oo 

uj(x,y,z,[)  =  Ujj>(z)exp(ikmx)exp(ikyy)exp(-icot),  (2) 

m  =  — oo 

m  =  +°° 

Vj(x,y,z,t)=  22  vm(z)exp(ikmx)exp(ikyy)exp(-icot),  (3) 

m  =  -oo 

and 


Wj(x,y,z,t)  =  wm(z )  exp(ikmx)  exp(ikyy)  exp(-icot), 


(4) 


where 


km  —  kx  + 


271171 


(5) 


where  kx  is  wavenumber  with  respect  to  the  x-axis,  ky  is  the  wavenumber  with  respect  to  the  y-axis,  cd  is  frequency, 
t  is  time,  i  =  V^T  and  j  is  either  1,  which  corresponds  to  the  backing  plate  or  2,  which  corresponds  to  the  acoustic  coating. 
The  functions  U$(z),  V$(z),  and  W$(z)  are  given  by 


U%(z)  =  A%ikm  cos[oc®z]+B®ikm  sin[a ®z]  +  C®  ^  sin[/?®z]-D®  ^  cos [/?®z] 

Pm  Pm 


],2- 

bP+ol 

Pm  ^  R(j) 

rm. 


sin[j8®z]-F® 


bP  +  ol 

Pm  ^  p(j ) 
lJm. 


cos  [^®z], 


(6) 


V®  (z)  =  A®  iky  cos[a®z] + B®  iky  sin[a®z]-C® 


pO)  _j_ 

Pm  i"  o(D 
rm. 


sin[^z]+D^ 


p(j)  ,  Kn_ 
Pm  ^  p(j ) 

lJm. 


cos[£®z] 


_£®  %T  Sin^mZ]  +  Fm%^  C0S^mZ] 

Pm  Pm 


(7) 
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and 


W$(z)  =  ~Am«m  sin[a*z] +B®a®  cos[a®z]-C$iky  cos[0®z]-D®iky  sin[0®z] 
+E®ikm  cos[0®z]+F®ikm  sin[0®z], 

where  a$  is  the  modified  wavenumber  associated  with  the  dilatational  wave  and  is  expressed  as 


(8) 


O) 

where  is  the  dilatational  wavenumber  on  the  layer  (j)  and  is  equal  to  co/cf,  where  is  the  dilatational  wave  speed  on 
the  layer  (j);  is  the  modified  wavenumber  associated  with  the  shear  wave  and  is  expressed  as 

=  (io) 

where  kf  is  the  shear  wavenumber  on  the  layer  (J)  and  is  equal  to  oo/cf,  where  cf  is  the  shear  wave  speed  on  the  layer  (J). 
The  dilatational  wave  speed  is  given  by 


cd  = 


(12) 


and  the  shear  wave  speed  is  determined  by 


Q  — 


(13) 


If  the  model  of  the  layers  of  the  structure  contain  a  structural  loss  factor,  then  the  dilatational  and  shear  wave  speed  are 
complex  quantities  as  are  the  modified  wavenumbers  a  and  /?.  (Structural  loss  factors  are  typically  used  to  model  acoustic 
coatings  and  sometimes  used  to  model  metallic  hull  behavior.)  Note  that  the  positive  value  of  radical  is  chosen  in  Eqs.  (9) 
and  (10).  Choosing  the  negative  value  would  result  in  the  identical  solution  because  the  solid  layers  are  finite  in  the 
z-direction.  In  Eqs.  (6)-(8),  A$,  C$,  D £$,  and  F$  are  the  unknown  complex  wave  propagation  coefficients  whose 

solutions  we  seek. 

The  acoustic  pressure  in  the  fluid  medium  is  governed  by  the  three-dimensional  wave  equation  and  is  written  in 
Cartesian  coordinates  as 


e2Pa(x,y,z,t )  d2pa(x,y,z,t )  C^PaiX.y.Z.t)  1  d2pa(x,y,z,t)  _ 

8x2  +  8y2  +  8z2  cj  dt 2  ’  1  1 

where  pa(x,y,z£ )  is  the  pressure  and  cf  is  the  real-valued  compressional  wave  speed  in  the  fluid.  The  pressure  field  in  the 
fluid  consists  of  an  incoming  incident  pressure  wave  (system  excitation)  applied  to  the  structure  at  the  top  (z=c)  and  an 
outgoing  scattered  wave  that  propagates  in  the  positive  z-direction.  Due  to  the  periodicity  of  the  system,  the  outgoing 
pressure  wave  has  the  same  spatial  form  as  the  displacement  fields,  thus  the  solution  to  Eq.  (14)  is  written  as 

pa(x,y,z,t)  =  P/exp(-iy0z)  exp(ikxx)  exp(i/<yy)  exp(-icot) 

m  =  +oo 

+  Y)  exp(iy,„z)  exp(ik,„x)  exp(ikyy )  exp(-iojf),  (15) 

m  =  — oo 

where  P/  is  the  magnitude  of  the  applied  pressure  field,  Rm  are  unknown  wave  propagation  coefficients,  and  ym  is  the 
modified  wavenumber  associated  with  the  compressional  wave  in  the  fluid  and  is  given  by 

ym  =  \/(®/c/)2-k^-k|  =  ^kj-kl-k 2.  (16) 

The  expression  given  by  Eq.  (16)  is  either  purely  real  or  purely  imaginary.  The  positive  root  of  Eq.  (16)  is  associated  with 
the  outgoing  radiated  compressional  wave,  and  when  this  term  is  purely  imaginary,  the  sign  of  the  argument  of  the  series 
terms  in  Eq.  (15)  is  negative,  and  this  corresponds  to  a  spatially  decaying  wave  in  the  positive  z-direction.  The  relationship 
between  arrival  angles  on  the  plate  of  an  acoustic  wave  and  the  x-  and  y-wavenumbers  is  determined  by 

kx  =  (co/Cf)sm(6)  =  kfS\n(6)  (17) 


and 


ky  =  flj-k2  sin(f/>),  (18) 

where  6  is  the  arrival  angle  of  the  acoustic  wave  with  respect  to  the  x-axis  and  < j>  is  the  arrival  angle  with  respect  to  the 
y-axis.  Note  that  the  expression  under  the  radical  in  Eq.  (18)  is  always  positive,  and  that  the  positive  root  is  used  for  the 
analysis.  The  value  of  Eq.  (18),  however,  can  be  negative  based  on  the  sign  of  </>.  A  value  of  0=0  =  0  corresponds  to  a 
broadside  or  boresight  wave  exciting  the  system. 
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The  displacement  of  any  individual  rib  in  the  z-direction  is  modeled  using  the  Timoshenko  beam  equation,  which  is 


cld4sz(x,y,t) 

dy4 


84sz(x,y,t)  d2sz(x,y,t)  p2Id4sz(x,y,t) 
dy2dt 2  +  p  at 2  +  Ck  at4 


=  -fz(x,y,t)+ 


EI  a2/z(x,y,t) 

CAk  dy2 


pi  &fz{x,y,t) 
gak  at2  ’ 

(19) 


where  sz(x,y,t )  is  the  rib  displacement  in  the  z-direction,  /z(x,y,t)  is  the  force  per  unit  length  in  the  z-direction,  E  is  Young’s 
modulus,  I  is  the  second  moment  of  area  about  the  x-axis,  p  is  the  density,  G  is  the  shear  modulus,  A  is  the  area,  and  k  is  the 
shear  coefficient.  Noting  that  the  rib  displacement  in  the  z-direction  is  equal  to  the  plate  displacement  at  z=a  in  the 
z-direction,  i.e. 


m  =  +  oo 

sz(x,y,f)  =  Y  win 1  <  a  )exp(  i  km x)  exp(ikyy)  exp(-  ia>  t), 

m  —  -oo 

yields  the  solution  to  the  z-directional  force  in  Eq.  (19)  as 

m  =  +oo 

fz(x,y,t)  =  Kz  Y  wm  ><a>  exp(ik„,x)  exp(ikyy)  exp(-ia>t) 

m  =  — oo 

with 


EIk^-(pI-\-(pIE/GK))co2k2 -Apco2  +  ( p2I/Gk)cd 4 
-1-(EI/GAK)k$  +  (pI/GAK)co2 


(20) 


(21) 


(22) 


where  I<z  has  units  of  stiffness  per  unit  length.  Note  that  the  case  when  the  denominator  of  Eq.  (22)  equals  zero 
corresponds  to  the  maximum  frequency  that  the  Timoshenko  beam  equation  is  considered  valid.  This  frequency  limitation 
is  discussed  in  more  detail  in  Section  5. 

The  rotation  of  any  individual  rib  is  modeled  using  the  torsional  wave  equation,  written  with  the  independent  variable 
as  the  angle  of  the  rib  as 


d2e(x,y,t)  d26(x,y,t) 

pJ^F2 - GJ  W2 


-my(x,y,t), 


(23) 


where  0(x,y,t )  is  the  clockwise  angle  of  rotation  of  the  rib  about  its  centroid,  my{x,y,t)  is  the  torque  per  unit  length  in  the 
y-direction  at  the  top  of  the  rib,  and  J  is  the  polar  moment  of  inertia  about  the  y-axis.  Rewriting  the  torque  expression  as 

my(x,y,t)  =  Q)/Z(x,y,t),  (24) 


where  h  is  the  height  of  the  rib  and  /X(x,y,t)  is  the  force  per  unit  length  in  the  x-direction,  and  assuming  that  the  rib  is 
always  at  a  right  angle  to  the  bottom  of  the  backing  plate,  gives  the  relationship 

0(x,y,t)  =  =  Y  ikmw£\a)  exp(ikmx)  exp(ikyy)  exp(-icot)  (25) 

m  =  -oo 


which  yields 


m  =  +oo 

fx(x,y,t)  =  I<x  Y  ifemW<J)(a)exp(ikmx)exp(ifcyy)exp(-i®t) 

m  =  -oo 


(26) 


with 


kx=f[w2p-k2Gl  (27) 

where  I<x  has  units  of  stiffness  and  the  caret  has  been  added  to  this  constant  because  it  has  different  units  than  the  stiffness 
constant  in  the  z-direction. 

The  displacement  of  any  individual  rib  in  the  y-direction  is  modeled  using  the  longitudinal  wave  equation  of  a  bar, 
written  as 


Ap 


dt 2 


dy2 


(28) 


where  sy(x,y,z)  is  the  rib  displacement  in  the  y-direction  and  fy(x,y,t)  is  the  force  per  unit  length  in  the  y-direction.  The  rib 
displacement  in  the  y-direction  is  equal  to  the  plate  displacement  at  z=a  in  the  y-direction,  thus 


sy(x,y,t)  =  Y  '/m1)(a)exP(ifcmX)exp(ikyy)exp(-i®t) 


(29) 
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and  this  yields  the  solution  to  the  y-directional  force  in  Eq.  (28)  as 

m  =  +oo 

fy(x,y,t)  =  Ky  V^\a) exp(ikmx)  exp(i/cyy)  exp(-icof)  (30) 

m  =  -oo 


with 


Ky  =Apco2  -AEky,  (31) 

where  I<y  has  units  of  stiffness  per  length.  Note  that  the  dynamic  beam  stiffness  expressions  Kx ,  I<y,  and  I<z  can  change  signs 
based  on  the  parameters  of  the  beam  and  that  they  can  be  complex  quantities  if  the  material  properties  of  the  beam  are 
complex.  This  typically  results  when  a  structural  damping  model  is  used  which  produces  a  complex  Young  and  shear 
modulus  of  the  beam.  For  the  remainder  of  the  paper,  the  exponential  function  with  respect  to  time  is  suppressed  in  all  the 
equations. 

The  solutions  to  the  wave  propagation  coefficients  are  determined  by  formulating  the  problem  using  13  boundary  value 
equations  written  in  terms  of  the  plates’  displacements  and  corresponding  forcing  functions.  The  system  is  loaded  on  the 
bottom  by  the  ribs;  thus,  the  normal  stress  and  the  shear  stresses  at  this  location  (z=a)  are  written  using  a  force  balance 
between  the  forces  in  the  ribs  and  the  bottom  of  the  lower  plate.  Because  there  are  an  infinite  number  of  ribs  and  they  each 
act  as  a  point  force,  a  Dirac  comb  function  is  used  to  distribute  the  forces  of  the  ribs  discretely  in  the  x-direction. 
The  normal  stress  is  written  as 


Tzz(x,y,a,t)  =  2i 


dfa(x,y,a,t )  +  diq(x,y,a,t)  +  dwi(x,y,a,t) 


dx 


dy 


dz 


+  2  Ah 


dw^(x,y,a,t) 

dz 


=  fz(^y.t)S(x-nL) 

n  =  —oo 

and  the  tangential  stresses  are  written  as 

zzy(x,y,a,t)  =  fa 
and 

Tzx(x,y,a,t)  =  fa 


dv\  (x,y,a,t)  +  dwi(x,y,a,t)' 


dz 


dy 


dwi(x,y,a,t)  diq(x,y,a,t) 


dx 


dz 


=  Y2  fy(x,y,t)S(x-nL) 


=  Y)  fx(x,y,t)8(x-nL), 


(32) 


(33) 


(34) 


where  S(x-nL)  is  the  spatial  Dirac  delta  function,  and  the  subscript  1  corresponds  to  the  backing  plate. 

The  interface  between  the  backing  plate  and  the  acoustic  coating  is  modeled  using  continuity  of  displacements  and 
stresses.  Continuity  of  the  stress  fields  at  the  interface  (. z=b )  yields 


and 


X\ 


fa 


fa 


dfa(x,y,b,t)  +  dv\(x,y,b,t)  +  dw!(x,y,b,t)' 


dx 


dy 


dz 


+  2  fa 


dwi(x,y,b,t) 


=  2? 


du2(x,y,b,t)  dv2(x,y,b,t )  dw2(x,y,b,t) 


dx 


dy 


dz 


+  2  fa 


dz 

dw2(x,y,b,t) 


dz 


dvi(x,y,b,t)  +  dw!(x,y,b,t)' 


dz 


dy 


dwi(x,y,h,t)  dUi(x,y,b,t) 


dx 


dz 


=  fa 


-.fa 


dv2(x,y,b,t)  +  dw2(x,y,b,ty 


dz 


dy 


dw2(x,y,b,t)  du2(x,y,b,t)' 


dx 


dz 


(35) 

(36) 

(37) 


where  the  subscript  2  corresponds  to  the  acoustic  coating.  Continuity  of  the  displacement  fields  at  the  interface  yields 

fa  (x,y,b,t)  =  u2(x,y,b,t)f  (38) 


vi(x,y,b,t)  =  v2(x,y,b,t) 


(39) 


and 


w-L(x,y,b,t)  =  w2(x,y,b,t).  (40) 

The  top  of  the  acoustic  coating  is  in  contact  with  the  fluid,  and  the  normal  stress  at  this  location  (z=c)  is  written  using  a 
force  balance  between  the  pressure  in  the  fluid  and  the  coating  as 


Tzz(x,y,c,t)  =  22 


du2(x,y,c,t)  dv2(x,y,c,t)  dw2(x,y,c,t) 


dx 

=  -Pa(x,y,c,t), 


dy 


dz 


+  2  fa 


dw2(x,y,c,t) 

dz 


(41) 
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where  pa(x,y,c,t)  is  the  pressure  field  in  contact  with  the  top  of  the  plate.  Furthermore,  this  interface  satisfies  the  linear 
momentum  equation,  which  relates  the  normal  acceleration  of  the  plate  surface  to  the  spatial  gradient  of  the  pressure 
field  by 

„  82w2(x,y,c,t )  _  8pa(x,y,c,t)  fA^ 

Pf  dt2  ~  8z  '  K  ' 

where  pf  is  the  density  of  the  fluid.  The  tangential  stresses  on  the  top  of  the  acoustic  coating  are  modeled  as  free  boundary 
conditions  and  are  written  as 


and 


Tzy(x,y,c,t)  =  fi2 


dv2(x,y,c,t)  dw2(x,y,c,t )' 

az  +  dy 


=  0 


(43) 


Tzx(x,y,c,t)  =  p2 


~dw2(x,y,c,t) 

ax 


du2(x,y,c,t)' 

dz 


=  0. 


(44) 


3.  Analytical  solution 

The  functional  form  of  the  displacement  fields  in  Eqs.  (6)-(8)  are  now  inserted  into  the  boundary  value  equations  given 
by  Eqs.  (32)-(41)  and  (43)-(44).  Additionally,  the  pressure  field  in  Eq.  (15)  and  the  interface  equation  listed  as  Eq.  (42)  are 
utilized,  and  which  directly  results  in  the  following  equations: 

m  =  +oo  m  =  +oo 

k  Y  ifcmU!„)(a)exp(ikmx)+l1iky  Y  V'<|)(a) exp(ikmx) 


n  =  +oo  m  =  +oo 


171  =  +°° 

+  (21+2Ju1)  Y  — exp(ikmx)  =  Kz  Y  l  Y  W(V(a)exp(ikmx)]d(x-nL), 


(45) 


n  =  — oo  m  =  —oo 


n  =  +oo  m  =  +oo 


m  =  +°o  jw(1)/ n\  m  =  +oo 

Hi  Y  dz  exP(ifemX)+^iifcy  Y  WjJ>(a)exp(ikmx)  =  Ky  Y  [  Y  V'<J)(a)exp(ikmx)]<5(x-nL),  (46) 


m  =  -oo 


n  =  — oo  m  =  — oo 


n=+oom=+oo 


m=+0°  m  =  +OO  Jrr(l)/„x 

Hi  Y  ikmW^)(a)exp(ikmx)+p1  Y  — exp(i kmx)  =  Kx  Y  t  Y  ikmW^\a)exp(ikmx)]d(x-ni),  (47) 


n  =  — oo  m  =  — oo 


k  Y  ikmU$(b)exp(ikmx)+kiky  Y  V(m\b)exp(ikmx)+(k+2Hi)  Y 

m  =  -oo  m  =  -oo  m  =  -< 

m  =  +oo  m  =  +oo  m  =  + 

=  k  Y  ikmU%\b)exp(ikmx)+l2iky  Y  V^\b)exp(ikmx)+(X2  +  2p2)  Y 


m  =  -oo  dz 

m=  +00dWm\b) 


dz 


m  =  +°°  dv^(b)  m=+oo 

Hi  Y  df  exp(ik,„x)+/i1ikv  Y  W^(b)exp(ikmx) 


m  =  -oo 

m=+0°dV%\b) 


=  V2  Y 


m  =  -oo 
m  =  +oo 


dz 


exp(ikmx)+p2iky  Y  W%>(b)exp(ikmx), 


m  =  +00  m  =  +OO  JrrCO/L^ v 

/(,  ifemW‘,1)(b)exp(ikmx)+/i1  — 5^  exp(ikmx) 


m  =  -oo 
m  =  +oo 


dz 


m  =  -oo 

m=+0°dU%\b ) 


=  Hi  ^  ikmW®(b)exp(ikmx)+jU2  ^  exp(ikmx), 

m  =  -oo  m  =  -oo 

m  =  +oo  m  =  +oo 

X]  Um\b)  exp(ikmx)  =  X  Ug\b)  exp(ikmx), 

m  =  -oo  m  =  -oo 

m  =  +oo  m  =  +oo 

Y  V^(b)exp(ikmx)=  Y  ^>(b)exp(ikmx), 

m  =  -oo  m  =  -oo 

m  =  +oo  m  =  +oo 

X  W<|)(b)exp(ikmx)=  X  W%\b)  exp(ikmx), 


exp(i/cmx) 

exp(ikmx),  (48) 

(49) 

(50) 

(51) 

(52) 

(53) 


m  =  -oo 


m  =  -oo 
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m  =  -oo 


and 


m  =  -t-oo  m  =  +oo  m  =  +  c 

22  E  ikmU®(c)exp(ikmx)  +  A2i/cy  ^  V<,2)(c)exp(ikmx)+(/l2+2yU2)  ^ 

m  =  -oo  m  =  -c 

m  =  +oo  / (jq2 p  \ 

+  E  ;  wm2)(c)  exp(ikmx)  =  -2P,  exp(ikxx), 

m  =  -oo  \  l'm  J 

m  =  +  oo  jw( 2)  ✓  x  m  =  +  oo 

P2  E  dV  exp(ikmx)+/(2iky  ^  W®(c)  exp(ikmx)  =  0, 


m=+oo  m=+oo  jf/(2)/  x 

p2  E  ikmW^)(c)exp(ikmx)+/i2  ^  exp(ikmx)  =  0. 


3dWjh2)(c) 

dz 


exp(ikmx) 


(54) 

(55) 

(56) 


m  =  -oo  m  =  -oo 

To  eliminate  the  Dirac  delta  function  present  in  Eqs.  (45)-(47),  the  Fourier  series  of  the  Dirac  comb  function  is 
written  as 


ln^°°  (\2nnx\ 

E  S(x~nL>  =  i  E  exp(— l 


■ 


Furthermore,  the  identities 


n  =  +oo 

m  =  +oo 

[  vC(a)  1 

E 

E 

V<!>(a) 

|  exp(ikmx) 

n  =  -oo 

m  =  -oo 

[  ikmw£\a)  J 

exp 


^i27inx^ 

vw! 


wl?\a)  V 

E  1  ^(a) 

i  knWl?\a) 


(57) 


E  exp(ik,)ix)  (58) 


are  also  applied  to  Eqs.  (45)-(47)  to  separate  the  embedded  double  summation  into  two  multiplicative  summations. 
A  proof  of  Eq.  (58)  is  given  in  Appendix  A.  This  modified  version  of  Eqs.  (45)-(47),  with  Eqs.  (48)-(56)  are  all  multiplied  by 
exp(-i kpX)  and  integrated  over  [0 ,L].  Because  the  exponential  functions  are  orthogonal  on  this  interval,  the  equations 
decouple  into  sets  of  m-indexed  equations,  each  one  expressed  as 


M  i  km  Urn  (a)  +  (fl)  +  (^1  +  2/r1 ) 


dwti\a)  Kzn 


dz 


=  f  E 


(59) 


P,  +P, ikyW^Ha)  =  (60) 

L  n  =  -oo 

H,ikmw£\a)+n,  =  ^Y^iknWPia),  (61) 

uz  L  n  =  -oo 

mmU$,\b)+MikyV£\b)+(h  +2/i0  dVV£>(b)  -22ikmUg>(b)-22ikyVg>(b)-a2  +  2fe)dW£(b)  =  0.  (62) 

P,  +P,  ikyw£\b)-n2  dV^ -fi2ikyW^(b)  =  0,  (63) 

P]  dU”f(b)  -H2ikmW%\b)-n2  dU"foh)  =0,  (64) 

Ug>(i»-U®(i»  =  0,  (65) 

V£\b)-VV>(b)  =  0,  (66) 

WS)(ft)-W®(ft)  =  0.  (67) 

22ikmU®(c)+22ikyV^(c)+a2+2(u2)^^  +  W®(c)  =  |  “2P''  (68) 

(69) 

and 

/i2ikmW®(c)  +  /i2^fe=0.  (70) 
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Note  that  the  left-hand  side  of  Eqs.  (59)-(70)  are  m  indexed  and  the  right-hand  side  of  Eqs.  (59)-(61)  are  n  summed. 
Physically,  the  left-hand  side  of  these  equations  model  the  acoustic  coating  and  backing  plate,  the  right-hand  side  of 
Eqs.  (59)-(61)  models  the  ribs,  and  the  right-hand  side  of  Eq.  (68)  models  the  incoming  acoustic  energy. 

The  functional  forms  of  the  displacements  given  in  Eqs.  (6)-(8)  are  inserted  into  Eqs.  (59)-(70)  and  the  resulting 
algebraic  matrix  equation  for  each  set  of  m-indexed  coefficients  is 

f  {p},  m  =  0, 

[A(km)]{xm}=  ^  [F(fc„)](x„)  +  <  Q  (71) 

where  [A (km)]  is  a  12  x  12  matrix  that  models  the  dynamics  of  the  backing  plate  and  the  acoustic  coating,  {xm}  is  the  12x1 
vector  of  unknown  wave  propagation  coefficients,  [F(/cn)]  is  the  12x12  matrix  that  models  the  dynamic  interaction  of  the 
ribs  and  the  backing  plate,  and  {p}  is  the  12x1  vector  that  models  the  incoming  acoustic  wave  acting  on  the  structure. 
The  entries  of  the  matrices  and  vectors  in  Eq.  (71 )  are  listed  in  Appendix  B.  Eq.  (71 )  is  now  written  for  all  values  of  the  index 
m  and  the  results  are  rewritten  in  global  matrix  form.  This  mathematical  process  is  previously  described  [20]  and  this 
results  in 


Ax  =  Fx  +  p, 


where  A  is  a  block  diagonal  matrix  and  is  written  as 


0 

0 


F  is  a  rank  deficient,  block  partitioned  matrix  and  is  equal  to 


A(/c_!)  0  0 


A(/<0)  0 
0  A  (/<!) 


F  = 


p  is  the  system  excitation  vector  and  is  written  as 


(72) 


(73) 


F  (k-i)  F(/<0)  F(/<i) 
F  (/<_!>  F(/<0)  F(  kA) 
F  (k-0  F(/<0)  F ‘(Id) 


p  =  [--  0T  pT  0T  •••f 


(74) 


(75) 


and  x  is  the  vector  that  contains  the  unknown  wave  propagation  coefficients  and  is  equal  to 

x  =  [  •  •  •  {x_i  }T  (Xo)T  {Xl  }T  (76) 

where 

(Xo)  =  {  Aq1)  Bf  C®  D®  E®  F®  A®  Bff  C®  D<2)  Cf  F®  }T.  (77) 

The  0  term  in  Eq.  (73)  is  a  12  x  12  matrix  with  all  zero  entries  and  the  0  term  in  Eq.  (75)  is  a  12  x  1  vector  with  all  zero 
entries.  The  solution  to  the  wave  propagation  coefficients  is  now  found  by  truncating  the  matrices  in  Eq.  (72)  to  a  finite 
number  of  terms  and  solving 

x  =  [A-F]-1p.  (78) 

Once  these  are  known,  the  displacement  field  of  the  system  in  the  spatial  domain  can  be  determined  using  Eqs.  (2)-(8). 
Furthermore,  the  stress  distribution  in  the  backing  plate,  acoustic  coating  and  the  scattered  acoustic  field  can  be  computed. 


4.  Model  validation 

The  elastic  model  that  has  been  developed  in  Sections  2  and  3  can  be  compared  to  a  fluid-loaded,  ribbed,  Bernoulli- 
Euler  thin  plate  model  that  has  been  previously  developed  [7-9].  This  will  provide  validation  of  the  model  for  low 
frequencies,  low  wavenumbers,  and  small  plate  thicknesses.  This  Bernoulli-Euler  model,  however,  only  incorporates 
flexural  wave  behavior,  making  the  model  assumptions  invalid  at  higher  frequencies  and  wavenumbers.  The  thin  plate 
model  has  one  degree-of-freedom  that  is  displacement  in  the  z-direction  and  is  a  constant  value  across  the  thickness  of  the 
plate.  This  displacement  equation  is  written  as 

m  =  +oo 

w(x,y)  =  Wmexp(ifcmx)exp(ifcyy) 

m  =  -oo 


(79) 
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and  this  displacement  field  can  be  determined  by 


where 


and 


In  Eq.  (82),  the  indexed  entry  is 


with 


{w}  =  [T+K+R]-1  {f}, 


{«/}  =  {•■•  W_,  Wo  W,  ...}T 


m= 


"-ft 


r_!  o  o 

o  To  0 

0  0  Ti 


1  1  1 
1  1  1 
1  1  1 


[R]=(^ 


Tn  =  +fy-pta>2  +(pfa>2 /iyn), 

Et3 


D  = 


12(1— u2)’ 


(80) 

(81) 

(82) 


(83) 


1 

i 

i 

o 

1 

o 

o 

o 

l<0kj  ■  ■  ■ 

(84) 

fclfc-1 

k  ifco 

ki  ki 

0  —  2Pj  0  • 

-}T 

(85) 

(86) 


(87) 


where  t  is  the  thickness  of  the  plate,  E  is  Young’s  modulus,  and  v  is  Poisson’s  ratio. 

Fig.  2  is  a  plot  of  the  magnitude  of  the  transfer  function  of  the  plate  normal  displacement  (w)  divided  by  the  amplitude 
of  the  applied  incident  pressure  field  (P/)  versus  spatial  position  in  the  x-direction  at  a  frequency  of  50  Hz,  an  x-direction 


s 


I 


Spatial  Position  -  ^-Direction  (m) 


Not  to  Scale 


I 


3> 


Fig.  2.  Transfer  function  of  z-direction  displacement  divided  by  incident  pressure  versus  spatial  location  in  the  x-direction  at  50  Hz  for  the  thick  plate 
system  ( _ )  and  the  thin  plate  system  ( • ). 
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Fig.  3.  Transfer  function  of  U^(a)/Pi  (•),  vH)(a)/Pi  (▼),  and  w£\a)/Pi  (■)  versus  index  number.  The  dashed  line  is  the  machine  precision  of  the 
computation. 


wavenumber  kx  of  0.0209  rad  m-1  (0=30°),  and  ay-direction  wavenumber  ky  of  0.0256 rad m-1  (0=45°).  The  solid  line  is 
the  elastic  plate  theory  developed  above  and  corresponds  to  Eq.  (4)  and  the  dot  symbols  are  the  Bernoulli-Euler  plate 
theory  and  corresponds  to  Eq.  (79).  This  example  was  generated  using  an  extremely  low  frequency  and  small 
wavenumbers  because  this  is  a  region  in  which  the  thick  plate  model  should  theoretically  agree  with  the  thin  plate  model. 
The  following  parameters  were  used  for  the  plate:  total  plate  thickness  t  is  0.01  m,  plate  densities  pi  and  p2  are 
1200  kgm-3,  Lame  constants  X\  and  X2  are  9.31  x  10s  Nm-2,  Lame  constants  pi  and  p2  are  1.03  x  108Nm-2,  fluid  density 
pf  is  1025  kgm-3,  fluid  compressional  wave  speed  cf  is  1500  ms-1  and  spatial  location  y  is  0.  The  following  parameters 
were  used  for  the  ribs:  width  is  0.0064  m,  height  is  0.1016  m,  Young’s  modulus  is  200  x  109Nm-2,  Poisson’s  ratio  is  0.3, 
density  is  7860 kg m-3,  shear  coefficient  is  0.833  and  rib  separation  distance  L  is  1.0m.  These  parameters  result  in 
computed  dynamic  rib  stiffnesses  Kx  equal  to  7.955  x  103 Nm-1,  I<y  equal  to  4.159  x  105 Nm-2,  and  Kz  equal  to 
5.005  x  105Nm-2.  For  the  thick  plate  model,  the  plate  interface  region  location  b  is  -0.006m  and  the  output  location  is 
z=-t/2,  which  is  -0.005  m.  These  parameters  produced  a  dimensionless  fluid  wavenumber  multiplied  by  plate  height 
( kjh )  equal  to  0.0021,  which  is  a  very  small  number  where  Bernoulli-Euler  plate  theory  is  considered  accurate.  The  elastic 
plate  model  was  truncated  to  15  modes  that  produced  a  120-by-120  system  matrix  and  the  thin  plate  model  was  truncated 
to  201  modes.  There  is  almost  total  agreement  between  the  two  models  over  all  of  the  abscissa  values.  To  illustrate  the 
convergence  of  the  individual  terms  from  the  series  solutions  in  Eqs.  (2)-(4),  each  term  is  plotted  versus  the  index  number 
in  Fig.  3.  In  this  plot,  the  square  markers  are  the  terms  W^\a)/PIt  the  circular  markers  are  the  terms  U^(a)/PIt  the 
triangular  markers  are  the  terms  V^}(a)/PIf  and  the  dashed  line  is  the  working  precision  of  the  computation  using  64-bit 
arithmetic  numbers.  For  this  validation  problem,  there  is  46.2  dB  of  falloff  for  [/^(a)/?/,  a  57.1  dB  of  falloff  for  V^}(a)/P/,  and 
a  57.8  dB  of  falloff  of  Wm}(a)/Pj  between  the  3rd  and  15th  mode.  For  comparison,  a  1/m2  convergent  series  has  a  28.0  dB 
falloff  between  the  3rd  and  15th  term.  Other  tests  of  these  series  using  various  parameters  resulted  in  similar  behavior, 
however,  as  frequency  increased,  more  terms  were  needed  to  insure  the  sum  of  the  series  has  converged. 


5.  An  example  problem 

An  example  problem  is  now  formulated  and  discussed.  This  problem  consists  of  a  0.0065-m  (0.25-inch)  thick  aluminum 
plate  coated  with  a  0.0254-m  (1-inch)  thick  urethane  polymer  coating.  This  example  was  generated  with  the  following 
plate  parameters:  aluminum  Lame  constant  is  5.11  x  10loNm-2,  Lame  constant  pi  is  2.63  x  10loNm-2,  density  pi  is 
2710kgm-3,  urethane  Lame  constant  X2  is  2.09  x  109Nm-2,  Lame  constant  p2  is  7.14  x  107Nm-2,  density  p2 
is  lllOkgm-3,  urethane  structural  loss  factor  is  0.05,  fluid-compressional  wave  speed  C/ is  1475ms-1  and  fluid  density 
pf  is  1025  kgm-3.  The  following  parameters  were  used  for  the  ribs:  width  is  0.0064  m,  height  is  0.1524  m,  Young’s  modulus 
is  70  x  109Nm-2,  Poisson’s  ratio  is  0.33,  density  is  2710 kgm-3,  shear  coefficient  is  0.833  and  rib  separation  distance  L  is 
1.0  m.  No  structural  loss  factors  are  used  for  the  hull  or  ribs  because  they  are  typically  metallic  for  a  marine  vehicle 
and  have  an  extremely  low  loss  factor,  however,  this  term  can  be  included  in  the  analysis  using  complex  material 
properties  for  these  components.  For  this  example,  the  material  properties  are  considered  constant  over  all  frequencies 
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that  the  analysis  is  undertaken.  Frequency  dependent  material  properties  can  easily  be  incorporated  into  the  model  if 
needed.  The  problem  is  investigated  from  three  technical  standpoints:  (1)  the  dynamic  response  of  the  structure,  (2)  the 
reception  of  an  acoustic  signal  from  within  the  structure,  and  (3)  the  scattered  pressure  field  of  an  acoustic  emission 
directed  at  the  system. 

Fig.  4  is  a  plot  of  the  normal  stress  in  the  z-direction  versus  x-spatial  position  and  frequency.  This  plot  is  at  the  midpoint 
of  the  combined  plate  system  which  is  z=  -0.0159  m  (the  top  of  the  plate  system  is  z=0),  the  value  of  y  is  0,  and  the 
incoming  acoustic  wave  has  wavenumbers  of  kx= 0  and  ky= 0.  The  scale  of  the  plot  is  power  in  decibels  and  the  units  are 
dimensionless.  In  addition  to  the  normal  stress  in  the  z-direction,  the  normal  stresses  in  the  x-  and  y-directions,  the  shear 
stresses  in  the  xy-,  xz-  and  yz-directions,  and  the  displacements  in  the  x-,  y-  and  z-directions  are  available  from  the  model. 
This  provides  a  three-dimensional  elastic  analysis  of  this  system.  Note  that  if  there  were  no  ribs  present  in  the  system,  this 
plot  would  be  constant  in  the  x-spatial  position  direction.  The  ribs  introduce  dynamic  effects  and  this  results  in  spatially 
varying  stress  and  displacement  fields.  This  amount  of  variation  is  dependent  on  the  stiffness  of  the  plate  system  with 
respect  to  the  stiffness  and  spacing  of  the  ribs. 

The  reception  of  an  acoustic  signal  in  the  urethane  can  be  accomplished  by  an  internal  array  of  sensors  embedded  in  the 
coating.  For  the  above  example  problem,  an  eight-element  linear  array  of  sensors  that  has  an  element-to-element  spacing 
of  0.2  m  and  the  first  sensor  located  at  x= 0.02  m  is  analyzed.  The  array  is  oriented  such  that  it  is  parallel  to  the  x-axis. 
The  length  of  this  array  makes  it  cross  one  rib  in  the  structure,  and  the  effect  on  the  response  of  this  crossing  must  be 
incorporated  into  the  analysis.  The  magnitude  of  the  response  of  two  points  at  a  distance  L  apart  are  identical;  however,  the 
phase  angle  must  be  adjusted  to  account  for  the  orientation  of  the  acoustic  wave  to  the  structure.  This  will  result  in  a  phase 
angle  difference  of  two  points  a  distance  L  apart  of  kxL  radians.  Fig.  5  is  a  plot  of  the  normalized  summed  array  response 
versus  arrival  angle  at  3600  FIz,  which  is  approximately  the  Nyquist  frequency  of  the  array.  The  summed  array  response 
was  computed  using 


N 

B(kx,co)  =  Y  77(x„,0, zs,co)  exp(-i ksxn\  (88) 

n  =  1 

where  77  is  the  field  variable  of  the  sensor  response,  N  is  the  number  of  sensors  in  the  array,  ks  is  the  steered  wavenumber 
of  the  array  with  respect  to  the  x-axis,  and  xn  is  the  location  of  the  nth  sensor.  In  Fig.  5,  the  solid  line  in  the  top  plot 
corresponds  to  an  array  of  z-direction  displacement  type  sensors  (typically  an  accelerometer)  and  the  solid  line  in  the 
bottom  plot  corresponds  to  an  array  of  z-direction  stress  type  sensors  (typically  a  hydrophone).  In  both  the  top  and  the 
bottom  plots,  the  dashed  line  corresponds  to  an  array  response  in  the  system  without  ribs.  In  this  example,  the  incoming 
acoustic  wave  has  wavenumbers  of  kx= 0  and  ky= 0.  The  plots  are  not  symmetric  about  the  steer  angle  of  0°  because  the 
array  is  not  symmetric  with  respect  to  the  structure.  For  this  example,  the  interaction  of  the  ribs  with  the  sensor  field  is 
significant  for  the  stress  sensors,  as  the  array  response  produces  higher  sidelobe  levels  than  an  ideal  array  would  produce. 
There  is  a  smaller  effect  of  the  ribs  on  the  summed  array  response  of  displacement  sensors,  although  there  is  still  some 
increase  in  the  sidelobe  levels.  Note  that  the  array  analysis  can  also  be  conducted  using  a  planar  array  embedded  in  the 
coating,  and  the  resulting  beam  pattern  would  be  two-dimensional. 


Fig.  4.  Transfer  function  of  z-direction  normal  stress  divided  by  incident  pressure  versus  spatial  position  in  the  x-direction  and  frequency  for  the  thick 
plate  system. 
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Fig.  5.  Summed  acoustic  response  versus  steer  angle  for  an  array  of  (a)  z-direction  displacement  sensors  ( _ )  and  (b)  z-direction  normal  stress  sensors 

( _ )  compared  to  thick  plate  system  without  ribs  ( _ ). 


In  underwater  structures,  the  scattered  pressure  field  is  an  important  quantity  for  quiet  operation.  The  scattered 
pressure  field  divided  by  the  incident  pressure  field  is  given  by  the  expression 

P  (x  V  0)1  m  =  +  oo  q 

5  'n'  -  =  exp(ikxx)exp(ikyy)+  JJ,  wm(0)  exp(ikmx)  exp(ikyy),  (89) 

11  m  =  -oo  l'm 


where  Ps{x,y,co )  is  the  scattered  pressure  field  in  the  spatial-frequency  domain,  the  first  term  on  the  right-hand  side  is  the 
reflected  pressure  field,  and  the  second  term  on  the  right-hand  side  is  the  radiated  pressure  field  caused  by  normal 
displacement  of  the  plate.  Fig.  6  is  a  plot  of  the  scattered  pressure  field  divided  by  the  incident  pressure  field  versus 
x-spatial  position  and  frequency.  The  incoming  acoustic  wave  has  wavenumbers  of  kx= 0  and  ky= 0.  The  scale  of  the  plot  is 
power  in  decibels  and  the  units  are  dimensionless.  The  scattered  pressure  field  also  has  spatial  variation  in  the  x-direction. 

The  above  simulations  (and  others)  were  run  using  all  three  nonzero  dynamic  stiffness  values.  These  were  compared  to 
simulations  where  I<x  =  I<y  =  0  and  this  resulted  in  nearly  identical  values  for  the  displacements,  stresses  and  scattered 
pressure  fields.  This  result  makes  sense  physically  because  the  dynamic  stiffnesses  I<x  and  I<y  act  on  the  in-plane 
displacement  fields  of  the  plates  while  I<z  acts  on  the  normal  displacement  field.  Because  the  plate  system  is  much  stiffer 
in-plane,  the  contribution  of  stiffeners  in  these  directions  will  be  minimal.  Utilizing  the  fact  that  the  dominant  beam 
stiffness  effects  are  in  the  z-direction  of  the  system,  an  upper  frequency  limit  of  the  model  can  be  determined.  For  the 
Timoshenko  beam,  the  upper  frequency  is  usually  given  by 


&>MAX  : 


For  a  rectangular  beam,  this  expression  can  be  approximated  by 

/max~^. 


(90) 


(91) 


where /max  has  units  of  FIz,  cs  is  the  shear  wave  speed  in  the  beam  and  h  is  the  height  of  the  beam.  For  the  example  problem 
presented  here,  this  upper  frequency  limit  is  computed  to  be  10,290  FIz.  The  plate  equations  without  the  ribs  have  no  upper 
frequency  bound  as  these  equations  of  motion  are  fully  elastic. 
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Fig.  6.  Transfer  function  of  scattered  pressure  field  divided  by  incident  pressure  versus  spatial  position  in  the  x-direction  and  frequency  for  the  thick  plate 
system. 

6.  Conclusions 

An  elastic  analytical  model  of  a  system  that  consists  of  a  ribbed  plate  covered  by  a  fluid-loaded  acoustical  coating  has 
been  derived.  This  model  was  developed  so  that  structural  acoustic  modeling  of  such  a  system  could  be  undertaken  at 
higher  wavenumbers  and  frequencies,  rather  than  the  previously  available  low  frequency  and  low  wavenumber  models. 
The  new  model  has  been  shown  to  agree  with  previously  developed  thin  plate  solutions  to  the  problem  at  low  frequency. 
An  example  problem  was  developed  to  illustrate  high-frequency  behavior  and  an  analysis  of  the  dynamic  response, 
acoustic  reception,  and  scattered  acoustic  field  is  included.  It  was  shown  that  for  this  example,  the  ribs  have  a  significant 
effect  on  the  high  frequency  dynamic  response  of  the  structure. 
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Appendix  A.  Double  summation  identities 


A  proof  of  the  double  summation  identities  used  in  Eq.  (58)  is  developed  here.  Starting  with  the  series  Gm  (where  Gn 
equals  either  V\4J}(a),  V^\a),  or  i/<m W^J)(a>),  the  following  expression  is  written  as: 


n  =  +oo  m  =  +oo 


Y  t  Y  Gm  exp(ikmx)] exp (—j—J  ■ 


n  =  — oo  m  =  — oo 


Using  the  definition  of  km  from  Eq.  (5)  yields 

n=  +  oo 


exp(ik*x)  Y 


'^4°°  _  ra  nmx\ 

Y  Gmexp(— — ) 


exp 


Expanding  the  m  indexed  series  results  in 

n  =  +oo 


exp(i/<xx)  ^ 


+  G_iexp 


i27i(-l)x 


+  G0  exp 


Zi27i0x\ 


^i27inx^ 


/i27ilx\ 


(v-^J+Giexp(v^J+---Jexp(v: 


Next,  expanding  the  n  indexed  series  yields 


exp(i/cxx)<  •  •  •  + 


i27r(— l)x 


Zi27i0x\ 


/i27ilx\ 


•  •  •  +  G_!  exp ^ - - - J  +  G0 exp —)  +  Gx  exp —)  +  • 


/i27inx\ 

~Trr 


exp 


l27T(-l)X 


(A.1) 


(A.2) 


(A.3) 
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\2n(-\)x 


f\2n0x\ 


f\2n\x\ 


+  G_,  exp {—j— )  +  Go  exp (  —  j  +  G,  exp j  +  ■  ■  ■  j  exp (  — 


i2n(-l)x 


( i2n0x\ 


+  G_!  exp ^ - - - J  +  G0 exp J  +  Gt  exp (  — - 


_  /i27ilx\ 


exp 


/i27iOx^ 

/i27ilx\ 

(— )+" 


Incorporating  the  outer  exponentials  into  the  inner  exponentials  gives 


exp(i/cxx) 


+  G_  i  exp 


i27i(-2)x 


+ 


+  G_i  exp 


L 

127T( — l)x 


+  G0  exp 


i27i(-l)x 


+  Gi  exp 


/i27i0x\ 


+  G_i  exp 


I 

/i27i0x\ 


/i27i0x\ 


+  G0 exp j— J  +  Ct  exp j- 


/i27ilx\ 


7i0x\  _  /i27ilx\  _  Zi27i2x\  1  1 

ir)+Coexp{—)+C'exp[—)+--  ]  +  ■"}• 


Summing  the  above  expression  diagonally  using  similar  exponentials  results  in 


exp(ikxx) 


+  c„ exp 


i27t(-l)x 


n  =  +oo 


n  =  — oo 


n  =  — oo 

Factoring  out  the  Gn  series  yields  the  expression 

n=^°°  r  /i27T(-l)x 

|  ■  ■  ■  -I-  caF 

n  =  — oo 

which  can  be  rewritten  as 


i27i(0)x\ 


+  £  GnexPpp  +  £  G„exPpp  + 


i27r(l)x\ 


exp(i/cxx)  Y  Gn 


V 


+  exp 


/i27i(0)x\ 


/i2rc(l)x\ 


^j+exp^J  +  ... 


n  =  +  oo  m  =  +  oo 


exp(ikxx)[  ]T  G„]  ]T  exp 


n  =  -oo  m  =  -oo  \  ~  J 

Finally,  incorporating  the  leading  exponential  into  the  m  indexed  series  gives  the  expression 

n  =  +  oo 


[  Y)  G"]  Y)  eXP (ikmX)- 


(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 


Appendix  B.  Matrix  and  vector  entries 

The  entries  of  the  matrixes  and  vectors  in  Eq.  (71 )  are  listed  below.  Without  loss  of  generality,  the  top  of  the  top  plate  is 
defined  as  z=c= 0.  For  the  [A kn]  matrix,  the  nonzero  entries  are 

au  =  (-2i[(a(n1))2+/c^+ky)]-2/r1(a(n1))2}cos(a}11)a), 
a  1,2  =  {-/Li[(ai,1))2  +  k2+k2)]-2JUi(a<11))2}sin(a<11)a), 
ali3  =  2i/i1^<„,)kysin(^1)a), 
a,,4  =  -2i  n^nky  cos 
a,, 5  =  -2i .n-itffkn  sint/l^a), 
ai.e  =  2i Hiftfkn  cos (jS^a), 
a2,i  =  -2i/r1a[11)/<y  sin(a^1}a), 
a2,2  =  2i^1a(n1)/<y  cos(a(n1}a), 
az.3  =  ~^[(P(n))2  +  k2„-kj]cos(P(Ja), 

a2,4  =  -tnKfiff+tf-fysinCffia), 

a2. 5  =  -2/(,kyk„  cos  (jS^a), 
a2,6  =  -2/q  kyk„  sin(/?'rl'Ja), 
a3,i  =  -2i/t1a(n,)k„  sinfa^a), 
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a3, 2  =  2iJu1aJJ1)k„  cos^'a), 
a3,3  =  2  jixknky  cos  (jS^a), 
a3,4  =2/ilk„/<ysin(/41,a), 

03,5  =  ^1  [(jSn1’)2— +  fcy]  COSC^^a), 

03,6  =  +fcy]  sin(jS<„1)a), 

a4,i  ={-li[(a<„1))2  +  k2  +  fe2)]_2jUl(af1))2}cos(a<„1)b), 
a4,2  =  {-^i[(an))2+^n+^y)]-2Ju1(a{11))2}  sinCa^b), 
o4,3  =  2iji^ky  sin  (j^b), 
a44  =  —  2  i /( ,  />  Jj1 1  ky  cos(  'b ), 

o4,5  =  -2i  Hxffikn  sin(//„"b), 

04,6  =  2i  Hxftn’kn  C0S(/^,"b), 
a4,7  =  {^2[(aS,2))2  +  k2  +  ky)]+2jU2(aJ,2))2)  cos(a(„2)b), 
a4,8  =  U2[(a®)2+k2+^)]+2/i2(a®)2}sin(a®b), 
o4l9  =  —  2iJu2/?(„2)ky  sin(/?(n2)b), 
o4,io  =  2i  n2P^ky  cos  (p^b), 

04.li  =  2ijU2^2)k„  sin  (jS^b), 

04, 12  =  -2i n2^kn  cos(jS<2)b), 
a5>  i  =  -2i^1a(n1)/<y  sin(a(n1}b), 
a5>2  =  2i^1a[11)/<y  cosCa^b), 

05,3  =  -jOiK^n’V+^-kylcos^^b), 

o5,4  =  -J4\ [(Pn^2  +  fc2  -fc2]  sin^^b), 

a5,5  =  -2/i, kyk„  cost fl'J’b), 

05. 6  =  -2  Hikykn  sim/^'b), 
a5>7  =  2i .fcoffiky  sin(aj,2)ib), 
a5  8  =  cos(aJj2)b), 

05,9  =  /i2[(iSn2))2  +  k2-k2]cos(/?!,2)b), 

05,10  =  ju2[0O2  +fen-fcy]  sin(jS!,2)b), 

05,11  =2  ti2kykn  cos  (ffib), 

05,12  =  2  Ji2kykn  sin(/?<2)b), 
a6,i  =  -2i/i,a<„,)knsin(a<11)b)> 
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a6i  2  =  2iiu1o41)/<n  cosiot^b), 
a6>3  =  2^,k„ky  cos(//„"b), 
a6A  =  2n-lknkySin($)b), 
a6.5  =  Mi  [(i?^)2-^  +  k2]  cos(p(fb), 
Oe.e  =  Mi  [(/J<n1))2-k2  +  k2]  sin  (jS^b), 
a6.7  =  2iJu2a<2)kn  sin(af,2jb), 
a6,s  =  -2i/z2aj,2)kn  cos(a(2)b), 
a6l9  =  -2  M2knky  cosijS^b), 
a6,io  =  -2  M2k„ky  sin(/?'„2>b), 

06,11  =  -M2KiSn2))2-fcn+fe2]cos(i?n2)b). 

06,12  =  -M2[0S!,2))2-k2  +k2]  sin  0?®b), 

o7,i  =  ik„  cosCa^b), 
a7  2  =  i/cn  sin(a^1}jb), 

“7,3  =  W  Sin(/?n'b)’ 

Pn 

07,4  =  COSO^b), 

Pn 

07,5=  ^n)+|r)j  Sin^b), 

oy.6  =  -  COS (#>), 

07,7  =  -ik„  cos(a},2)b), 
a7,8  =  — ikn  sin(a®b), 

07,9  =  sin(/?j,2)b), 

Pn 

ai.w=l^cos(p(n2>b), 

Pn 

07.il  =  -  ^!,2)+  sin (^2)b), 

a7.i2  =  (Pn}+-^j  c°s (/?<2)b), 

a8,i  =ikycos(aj!1)b), 
a8,2  =  iky  sin(a'r'Jb), 
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a8.3  =  -  {Pn’+p-^j  sin (#>), 

asA=(ff+^Jcos(Wb), 

a8,5  =  sin ($>), 

Pn 

as,6=||cos0Cf>), 

Pn 

as,  i  =  -iky  cos(a*,2)b), 
a8i8  =  -iky  sin(a<2)b), 

as.9=  ^n2)+^)j  sin(/?(n2)b), 

as,  10  =  -  ^2)  +  ^  cos (ffi’b), 

as.ii  =  ^  sin(^2)b), 

Pn 

as,  12  =  cos(/?'„2’b), 

Pn 

a9,  i  =  -a™  sinCa^b), 
ag,2  =  a!,1’  cosfa^b), 
a9,3  =  -iky  cos(^fb), 
a9A  =  -ikysin(/?;’'b), 
a9,5  =  ikn  cos  (/5<„1)b), 
a9,6  =  ik„  sim/J^'b), 
a9,7  =  a(2)  sin(aj,2)b), 
a9,8  =  -a(2)  cos(aj,2)b), 

a9,9  =  iky  cos(/?'„2lb), 

09,10  =  iky  sin(/42)b), 
a9,n  =-ik„  cos(j3<2)b), 
a9,i2  =  -ik„  sin($,2)b), 

a10,7  =  -  ^2  [(«n2)  )2 + k„ + k2 )]  -  2/x2  (a(2))2 , 

C02P/«n2) 

iyn 

-co2pfky 


a10,9  = 


7, 
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aio.io  =  -2iyU2/?„2)ky, 

co2Pfkn 

aw'u  =  ^T’ 

a\0,U  =  24*2 

an, 8  =2\fi2c^)ky, 
all,9  =  —^2  (/42)  +  kn  — ky)» 
an, 11  =  —  2/i2/<n^y* 

ai2,s  =  2i/i2aJj2)/cn, 

ai2,9  —  2/^2  ^n^y* 

and 

Ql2.11  =A»2[(A2))2-kn+fcy]' 

The  nonzero  entries  of  the  [F(/cn)]  matrix  are 

/i,i  =  -(Kz/L)aJ,1)  sinfa^a), 

/i,2  =  (/<z/t)a(„1)  cos(ccl„1,a), 

/i,3  =  -(I<z/L)iky  cos(/?(„1)a), 

/i,4  =  -0<z/L)iky  siiK/^'a), 

/1.5  =  (Kz/L)ik„  cos(/?'„"a), 

/i,6  =  (I<z/L)ikn  sim/^'a). 

/2.1  =  (Ky/Lnky  costaj^a), 

/2,2  =  (Ky/L)iky  sinCa^a), 

/2,3  =  -a<y/i)  +  pjj  sintffa), 

h.4  =  (Ky/L)  cos(/f<1)a), 

f2,5  =  -( Ky/L sintf^a). 

Pn 

/2,6  =  (Ky/L)^cos(/f<1)a), 

/3.1  =  -(iknKx/L)a(^  sinloc^a), 

/3,2  =  (ikn^x/LJa^  cosia^a), 

/3.3  =  -(ik„Kx/L)iky  cos(//„"a), 

/3,4  =  ~(ik„Kx/L)iky  sin  (jS^a), 

/3,5  =  (iknKx/L)ikn  cos  (fifa), 

and 

/3.6  =  (iknKx/L)ik„  sin(/?'rl'Ja). 
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The  entries  of  the  p  vector  are 

p  =  [0  00000000  —2 Pi  0  0]T 
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